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Abstract. For more than three years now we have been conducting a spectroscopic 
survey of detached eclipsing binaries (DEBs) from the All-Sky Automated Survey 
(ASAS) database. Thousands of high-resolution spectra of >300 systems were secured, 
and used for radial velocity measurements and spectral analysis. In our sample we 
found a zoo of multiple systems, such as spectroscopic triples and quadruples, visual 
binaries with eclipsing components, and circumbinary low-mass companions, including 
sub-stellar-mass candidates. 


1. Introduction 

Due to the possibility of obtaining absolute values of a number of fundamental stellar 
parameters, detached eclipsing binaries (DEBs) are one of the most important objects 
in astronomy. However, having them in triple or higher order stellar systems is even 
more exciting. In addition to the many astrophysical applications of the “lone” debs, 
like testing the stellar structure and evolution models or building relations between 
stellar parameters, multiples can be used for testing the formation theories, dynamics, 
general relativity, and they also provide more stringent tests for the evolution models 
than binary stars. 

In our ongoing observational program that targets eclipsing binary stars from the 
ASAS survey (Pojmariski 2002) we have secured about 3000 high-resolution optical 
spectra of roughly 300 systems. We use this data to measure radial velocities (RV) 
of the binaries, to later combine them with photometry obtain a full physical model 
of a given object. In this sample we have identified a number of triple and quadruple 
systems, that contain an eclipsing pair. Of all the ~300 objects from our program, the 
multiples constitute 22 ± 2 % of the systems. We can divide them into two groups: 

• Systems with visual, resolved companions. 

• Systems with companions detected spectroscopically, which can produce: 

- RV trend of the eclipsing pair 

- the third set of spectral lines. 
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2. Eclipsing binaries with visual companions 
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Figure 1. Examples of visual companions to eclipsing binaries. Left: ASAS-011 
from Helminiak et al. (2012). Separation is 1.4 asec. Right: ASAS-05 (unpublished) 
observed with the AO 188 system at the Subaru telescope. The DEB is the star to the 
right. Separation is 0.385 asec, and the position angle is 62° (image rotated). 

These are systems where two stars are resolved (sometimes only with the adap¬ 
tive optics), and one of them shows eclipses. The first such system we identified was 
ASAS J011328-3821.1 (ASAS-011; Helminiak et al. 2012), another example is ASAS- 
05, (Fig. 1) which will be presented later. We take spectra of both stars - the DEB 
and the companion - and measure the RVs. We perform the standard steps of DEB 
modeling, obtaining the full set of parameters, i.e. masses, radii, temperatures and lu¬ 
minosities. The photometric information (brightness, colors) is also available for the 
third star. Additionally, if the companion is a single star and its RV suggest it’s grav¬ 
itationally bound to the DEB, we can perform a spectral analysis (for example with 
SME; Valenti & Piskunov 1996), from which we get the values of the companion’s 
T e ff , logg, and [M/H]. It is reasonable to assume that the binary itself has the same 
metallicity. This is important during isochrone fitting, to solve the age-metallicity de¬ 
generation, and better constrain the age and evolutionary status of the system. We thus 
have information about the masses and radii of two stars, but luminosities, temperatures 
and gravities of three, which constrains the age even more. 


3. Hierarchical triples - RV trends and the third light 

Most of the multiples we found show a trend in their radial velocities, which is a mani¬ 
festation of their motion around the center of mass, common with the outer companion. 
We thus deal with hierarchical triples, majority of which have the tertiaries fainter than 
the components of the DEB, and not seen in the spectra. Two such trends are presented 
in Fig. 2, for a case of a long- and short-term variations (the period of the outer orbit is 
known for the latter). 
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Figure 2. Examples of RV trends induced by the outer companion. The points are 
the values of the systemic velocity of the inner binary, calculated under the assump¬ 
tion of a constant mass ratio. Left: long-term case, when the period is unknown. 
Right: short-term case, outer orbit period is 168.1 d. 

The more interesting case is when we see the the third light in the spectra, we 
can trace velocity change, and it shows a motion around the common center of mass 
(Fig. 3). In such situation, we can treat the trend of the DEB and the movement of the 
tertiary as they were two velocities of a long-period SB2 1 . If the outer period is known, 
because it’s short enough or from other sources, we can perform the orbital fit, and 
obtain M 3 sin 3 (/ 3 ) and (Mi + Mf) sin 3 (/ 3 ). Even if the period is unknown, the mass ratio 
M 3 /(Mi + M 2 ) can be easily found, and since we know Mi and M 2 (from the DEB’s 
orbit), we get the third mass immediately. Combining M 3 with M 3 sin 3 (/ 3 ) gives us the 
inclination of the outer orbit, with respect to the plane of the sky. For the full geometry 
of the system, one needs the mutual inclination of the two orbits, which in principle can 
be estimated from the eclipse timing (Rappaport et al. 2013). Additionally, the light 
curve modeling gives us the fraction of the third light in the system, so the luminosity 
and colors of the tertiary, and from the colors we can estimate the temperature. 

In the case presented in Figure 3 we have the outer period (P 3 = 415.8 d) and 
eccentricity (^3 = 0.3) known from the eclipse timing. The inner period is 3.421 d, the 
masses measured are: Mi = 1.02, M 2 = 0.96 and M 3 = 0.80 M 0 , and the inclination 
is Z 3 = 68.5°, which means that the inner and outer orbit are non-coplanar. In fact, the 
mutual inclination, from the eclipse timing, is estimated to be 28°, so not enough to 
induce Kozai mechanizm (Kozai 1962). In most of our cases we can’t perform such 
an analysis, because we don’t see the third star in the optical spectra. For this reason, 
in 2014 we have started an infra-red spectroscopic campaign on the Subaru telescope. 
Having similar results for more systems will allow us, to study the angular momenta 
distribution in multiples, or to build a statistic on the mass ratios of the tertiaries to the 
inner binaries, which will be important for the formation theories. 


4. Double-double systems 

We have also noticed cases where the RVs of the third stars show a short (few or few 
tens of days) periodicity, or the DEB’s period (from the photometry) is in agreement 


Parameters related to the third star or the outer orbit have the index 3, and to the inner binary (DEB) - 
indexes 1, 2 or 12. 
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Figure 3. Example of a triple system where all stars can have their RVs measured. 
Left: the eclipsing binary (inner) orbit with the trend removed. Right: The trend of 
the binary (filled symbols), the RVs of the tertiary (open symbols) and the orbital 
fit to these measurements. The outer orbit’s period and eccentricity were known 
beforehand and held fixed. 


with the motion of only one component in the spectrum, and the other two sets of lines 
show a different period. We thus deal with two spectroscopic binaries, either eclipsing 
SB2 + SB1, or eclipsing SB1 + SB2. In fact, the two systems from Fig. 1 are such 
double-doubles, ASAS-05 beeing the case of SB1 + SB2. Again, to spot the unseen 
components we need to observe in the infra-red. However, below we present the only 
case from our sample so far of an SB2 + SB2, i.e. a system where four stars have their 
RVs measured from optical spectra. 
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Non-eclipsing binary orbit, trend removed 
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Figure 4. Radial velocities of the ASAS-81 system. The eclipsing binary orbit 
(top left), its long-term trend (top right), the non-eclipsing binary orbit (bottom left), 
and its long-term trend (bottom right). Masses of all four stars and the non-eclipsing 
orbital inclination can be measured. 
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This system (ASAS-81, unpublished) is composed of 1.47-day eclipsing, and 7.18- 
day non-eclipsing pairs 2 , both showing measurable, opposite-sign RV trends, meaning 
their motion around the common center of mass (Fig. 4). The analysis of the 1.47- 
day DEB gives absolute masses of Mi = 1.19 and M 2 = 1.04 M 0 , and the orbital 
solution for the non-eclipsing pair gives M 3 sin 3 (/ 34 ) = 0.21 and M 4 sin 3 (/ 34 ) = 0.16 M 0 
(#34 = 0.781) and a small eccentricity ~0.005. The scales of the two trends imply 
(M 3 + M 4 )/(Mi + Mi) ^ 1.355, which later gives M 3 = 1.69 and M 4 = 1.32 M 0 . With 
such masses, the non-eclipsing binary’s inclination is 134 = 30°, which means the two 
“small” orbits are non-coplanar. The second trend, which is much better measured than 
the one of DEB, reaches 1.2 km/s in less than a year. One can thus suspect, that the 
“long” orbital period is of the order of single years, maybe ~10 yr, and over 5 solar 
masses are packed within < 10 AU. We are thus dealing with a system that is very 
interesting from the dynamical and point of view, and important for the multiple stars 
formation mechanisms. 

As stated before, most of the double-double systems we found contain a cold com¬ 
ponent, not seen in the optical spectra, which is also the case for the aforementioned 
ASAS-05. However, observations done at the Subaru telescope with the Infra-Red 
Camera and Spectrograph (IRCS) revealed a probable orbital solution (Fig. 5). If con¬ 
firmed, the IRCS measurements would imply qyi - 0.54, and masses 1.13 and 0.61 M 0 . 
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Figure 5. Radial velocities of the ASAS-05 system. Left: The eclipsing SB1 orbit 
with initial results for the secondary (blue). Right: The non-eclipsing SB2. 
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